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Abstract 

Low-energy thermal eqnilibrinm states of strongly conpled J\f = 4 snpersymmetric 
Yang-Mills (SYM) theory on a three-sphere are nnstable with respect to flnctnations 
breaking the global 50(6) R-symmetry. Using the gange theory/gravity correspon¬ 
dence, a large class of initial conditions in the R-symmetry singlet sector of the theory 
was been identified that fail to thermalize [1,2]. A toy model realization of such states 
is provided by boson stars, a stationary gravitational configurations supported by a 
complex scalar field in AdAs-gravity. Motivated by the SYM example, we extend the 
boson star toy model to include the global 50(6) R-symmetry. We show that sufficient 
light boson stars in the R-symmetry singlet sector are stable with respect to linearized 
fluctuations. As the mass of the boson star increases, they do suffer tachyonic insta¬ 
bility associated with their localization on 5®. This is opposite to the behaviour of 
small black holes (dual to equilibrium states of W = 4 SYM) in global AdS^: the latter 
develop tachyonic instability as they become sufficiently light. Based on analogy with 
light boson stars, we expect that the R-symmetry singlet nonthermal states in strongly 
coupled gauge theories, represented by the quasiperiodic solutions of [2], are stable 
with respect to linearized fluctuations breaking the R-symmetry. 
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1 Introduction and summary 


Consider maximally supersymmetric A/” = 4 SU {N) supersymmetric Yang-Mills theory 
on a three-sphere in the planar limit (iV —)■ oo, Qym 0 with Qym^ constant) 
and at large’t Hooft coupling, Qym^ 3> 1. In this limit the theory is best described by 
its holographic dual [3,4] — type IIB supergravity in asymptotically AdS^ x space- 
time. The global symmetry of the flve-sphere geometrizes the R-symmetry of the 
SYM. To simplify the discussion, we focus on S'^-invariant initial conflgurations of the 
SYM, and their evolution, consistently described within supergravity approximation. 
The vacuum of the theory has a Casimir energy 



where L is the radius of the S^. An initial state with the energy 


B = (1 + £) 


E ^ E 

^vacuum ^vacuum ? 


if it equilibrates in the future, is described by a Schwarzschild black hole in AdS's xS\ 
with the entropy S, the temperature T, and the size^ r+ given by 



(1.3) 


^The size is defined as a radius of the measured at the black hole horizon. 
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The dynamical equilibration of a large class of such SYM states, with an additional 
assumption that the dynamics occurs in the R-symmetry singlet sector^, has been 
discussed in [5,6]. Once a black hole becomes sufficiently small (light), 

^ < 0.440234 ^ e < CeHi ^ 2.53616, (1.4) 

it suffers from the Gregory-Laffamme (GL) instability [7] towards localization on 
[8,9,11]. The latter implies, in particular, that any initial condition (off the equilibrium) 
of the SYM with the energy less than ecru can not equilibrate within the R-symmetry 
singlet sector; in other words, the R-symmetry must be spontaneously broken in the 
approach to thermal equilibrium for any sufficiently low-energy state. First attempts 
to construct thermally equilibrium states in A/” = 4 SYM with broken R-symmetry 
were undertaken in [9]. 

Here, we point to another possibility regarding the low-energy dynamics of A/” = 4 
SYM. In [5] (and further extended in [10]) it was pointed out that evolutions starting 
with initial data close to a single mode in AdS did not collapse. The class of non¬ 
collapsing solutions was extended in [1,2]. A typical representative of this class is a 
boson-star [1] or a boson-star-like [2] (quasi-periodic^) configuration that is character¬ 
ized by a broad distribution of a scalar prohle in AdS and the dynamical evolution 
such that the nonlinear dispersion of the scalar energy-density overcomes the focusing 
effects of gravity in the AdS-cavity. Since the dynamics of [1,2] was discussed entirely 
in AdSd+i-, it necessarily occurs in the R-symmetry singlet sector of the holographically 
dual gauge theory. It is natural to expect that unlike evolutionary trajectories that end 
up in the thermal state, the dynamics of the nonthermal states of [1,2,5,10] is consis¬ 
tently restricted to the R-symmetry singlet sector — these states are stable with respect 
to localization in the compact manifold of the full ten-dimensional gravitational dual. 
The reason being the widely distributed energy-density prohle of the scalar helds sup¬ 
porting the solution that shuts-off the tachyonic instabilities observed for small black 
holes in AdS^ x S^. In this paper we present some evidence that such a scenario is 
indeed realized. Thus, we are led to conjecture that there is a large class of nonthermal 
low-energy states in strongly coupled gauge theories with unbroken R-symmetry. 

The rest of the paper is organized as follows. In the next section we describe a toy 
model of asymptotically AdS^ x holographic correspondence which supports boson 

^The full symmetry of the dynamics is thus S'0(4) x S'0(6). 

recent paper [12] argues that these solutions are anchors of the AdS-stability islands. 
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stars. We construct boson star solutions first in the effective five-dimensional descrip¬ 
tion, and further generalize the model to a ten-dimensional setting, where the compact 
manifold is a five-sphere and the effective five-dimensional negative cosmological con¬ 
stant is produced by the self-dual five-form flux. In section 3 we study stability of the 
ground state boson stars with respect to the linearized fluctuations breaking the SO{6) 
symmetry of the five-sphere. We find that boson stars are indeed free from the tachy- 
onic instabilities, provided they are light enough. We emphasize that this is opposite 
to the fate of smeared small black holes in AdS^ — we expect that all states in A/” = 4 
SYM with vanishingly small energy and unbroken /^-symmetry are non-equilibrium. 
We conclude in section 4. 

2 Boson stars in a holographic toy model 

Boson stars are stationary gravitational solutions supported by a complex scalar field 
stress-energy tensor. In asymptotically AdS space-times they were originally discussed 
in [13]. We begin with boson stars in asymptotically AdS^ space-time supported by a 
massless complex scalar field, and then extend the model to a ten-dimensional setting. 

2.1 Five-dimensional perspective 

Consider an effective action"^ 

(^5 + 12 - 3 ( 90 ( 90 ) , ( 2 . 1 ) 

IottGat 9_v(5 

where Gn is a five-dimensional Newton’s constant, 0 = 0i -f i02 is a complex scalar 
field, and 

M5 = dM5XI, dM5 = RtxS\ X= {y e [0,1]} . (2.2) 

Adopting the line element as 

dsl = - f-ae-^^dO + —— + (1 - y)dn^^ , (2.3) 

y \ 4|/(1 -|/)a J 

where dSL^ is the metric of unit radius and a{y) and diy) are scalar functions of 
the radial coordinate y describing the metric, and further assuming that the complex 
scalar field varies harmonically 

4>i{yA) + i4>2{yA)=p{y)0‘^\ (2.4) 

"‘We set the curvature radius of the asymptotically AdS^ solution to L = 1. 
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the equations of motion describing the boson star form the following system of ODEs: 


0 =p" + 


{2y-l)a-y + 2 , 




1 /( 1 /- l)a 


V 


iv(y-l)a^ ’’’ 


g2(5^2^2 

0=6' + 2y{y-l){pf- 


0=a' + 2y{y-l){p'fa + 


2a^ 

{2 — y){a — l) 


yiy - 1) 


2a 


A physically relevant solutions to (2.5) must satisfy: 
■ asymptotically at the AdS boundary, he., ?/ —)■ 0 +: 


P =Po y^ + 0 {y ^), 

a =1 + 02 + 02 y^ + ^{y ‘^)) 

d = = 2po y^ + Oiy^); 


(2.5) 


( 2 . 6 ) 


■ at the origin of AdS, i.e., z = 1 — y ^ 0+: 

p =Po - z + 0 ( 2 ^), 

o=i-i(<i;)>;)V 2 +o(j=), (2.7) 

^ = 111 4 - ^(dsnpS)v z + o(z^-). 

We compute the mass M (x E — E^acuum of the boson star as 

M= f dy^-j^ (4|/(1-|/)a^(p')^+ cu^e^V) , (2-8) 

Jo y 

and its charge Q 


Q= dS^ / dyy/^J ^, = ig^'' , 

Js^ Jo 


'53 

^2 


n -2 I ^ il-y)e^P^UJ 

Q =271 / dy - 5 - 

Jo y a 


(2.9) 


Note that given po, or alternatively the charge Q, the numerical boson star solution is 
determined by 4 parameters. 


{ca , 02 ) Pq j dg } , 


( 2 . 10 ) 


which is precisely the order of the system of ODEs (2.5). 
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Figure 1: Frequency uj and mass M of j = {0,1,2} (blue,red,green) boson stars as 
a function of charge Q. The dashed lines represent perturbative relation between the 
mass and the charge (2.11). 


The spectrum of light boson stars, i.e., for pq -C 1, is given by 


Pj = Po ‘^^- 2^1 (-i ) 4 + j ; 2 ; I - y) + 0{pl) , = 4 + 2j + 0{pl) , 

J + ^ 

Oj = 1 + 0{pI) , 5j = 0{pQf , 

87r2p2 


Qj ~ (■ 


(j + 2)2((j +2)2-1) 


+ o(p‘), M, = Q, + o(q 5) = ^4 % + o(Q?). 


TT^ 


TT^ 


( 2 . 11 ) 


where j = 0,1, • • • is the index. For general po (or Q) the boson stars can be found 
numerically, solving (2.5)-(2.7) with a shooting method as developed in [14]. The 
results of this analysis are presented in hgure 1. 


2.2 Ten-dimensional model 

We would like to generalize the effective action (2.1) to a ten-dimensional setting where 
we would be able to study the stability of the boson star solutions of (2.1) with respect 
to linearized fluctuations breaking the global symmetry of the compact manifold (we 
choose the latter to be S^). Our starting point is type IIB supergravity, where only 
the metric and the Ramond-Ramond hve-form 7^(5) are turned on. In this case 
the equations of motion take the form: 

G™ ^ = 0, <;F(„ = o, f,„=*f,5,. (2.12) 
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We now add a complex scalar field 0 = 0i + i(j) 2 , modifying (2.12) as follows 

C"’ ^ 41“’ - = 0. ^ 

(iF(5) = 0 , F(5) = *F(5) . 

We use (2.13) to model ten-dimensional boson stars. The latter equations of motion 
are obtained in a toy model, where type IIB supergravity Lagrangian is supplemented 
with a complex scalar field: 


CiiB(x. [ -) ^ V-fi'(io) { ^(lo) - 3(90(90H- 

JMw \ / JMw \ 

(2.14) 

We are interested in the most general ansatz describing solutions with 5'0(4) x S'0(5) 
isometry®. To obtain an explicit expression for the equations determining such solutions 
we first fix the reparametrization invariance such that 



Qte — 9x9 — 0 . 


(2.15) 


We can thus write the line element as®, 

ds^Q = —cl {dtY + ‘2^gtx dtdx + c\ [dxY + Cg {dVL^Y + c\ {dOY + Cg sin^ 6{dVLiY , 

F( 5 ) = (cKo dO + ai dt + a 2 dx) A dVl^ + (cks d9 A dt + d9 A dx + as dt A dx) A dVL^ , 
Ci = Ci(t,x,9), gtx = gtx{t,x,9), ai = ai(t,x,9), 

(2.16) 


where dQ^ is a volume form on a unit radius and dQ^ is a volume form on a unit 
radius Next we can eliminate {ao,ai,a 2 } by imposing 5-form self-duality. 


ao — 


C4C5a5 sin^ 9 

+ gl ’ 


0.1 — 


(a3C2 - a4,gtx)4 sin^ 9 


C4c|a/ 


+ 9i 


cfc: 


l'-2 


c4c|ycy| + gi 


(2.17) 


The resulting equations constitute a system of partial differential equations resulting 
from the eight non-trivial Einstein equations (the particular expressions are rather 
involved and we will not present their explicit form at this point): 


^(10) _ ^(10) _ ^(10) _ ^(10) _ ^(10) _ PI 

'^tt — '^xx — *-^0303 ~ ^99 ~ '-’'O4O4 ~ ^ ’ 


Gif = G)f = g3"> = 0 , (2.18) 


®We follow here the general discussion in [11]. 

^Expression for the RR five-form is the most general decomposition that preserves S'0(4) x SO{5) 
symmetry of the ansatz. 
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together with the five-form equations: 


0 =dea‘i + 4 cot( 6 ')Q ;3 - ^ d^. In ■ 


0:50405 


cfci + gtx Cl 


0 =(9004 + 4 cot ( 0)04 — 

C2 


clcl + 9tx 

(^59tx^4. 


I a /cjcl + 


44 + 91 


0^544 o 1 

44 + 9tx * 


dt In 


9tx 

O 5 C 4 C 5 


03 dg In 


C 1 C 4 C 3 

C 24 


4 




9tx 


4^/44^x 44 + 91^4 


,2 ’ 


C 15 C 4 C 5 


4444 +9u 


— aA dg In 


C4C| 


„2 2 

01 O 2 


H— TT~4~T^^ 

44 + 9 + 


Cl 


+ 


^b9tx^4 

44 + 9 I''" 


dr In 


O 5 C 4 C 5 


+ 


O 3 C 2 


df. 


9tx 


clAA+l + gtx A^cl + g"^^ 4 ’ (2.19) 


0 =( 9 t 03 - asdt In 


C44V44 + 9tx «4C? 


„2 4 
C 2 C 5 


dr In 


0 : 4+4 

C 2 C 3 C 4 


Osgtx o , ct3c| 
— ^dxln^^ 


Co 


C 3 C 4 


dx\n 


44 + 9 ‘tx 


+ 


0491x^1 


dr In 


9tx 


Cl 02 cl{clcl + g'^x) Cl 02 ci 


049 tx 1 <^405 

- dt In 

C 3 C 4 


dt In 


2 2 I 2 

cfc^ + gtx C 1 C 2 
0 =dga^ — dt0i4 + dxOx,. 


Boson star solutions of (2.18) and (2.19) with S0{6) symmetry take the form; 


03 — 04 — 0 , 05 — 


sm^ xe 

cos^ X 


^-5{x) 


Cl = 


a uc 


cos a; 


C2 = 


4a{x) cos X ’ 


gtx 0, 


C 3 = tanx , C 4 = C 5 = 1 , (j) = p(a;)e*“ , 


( 2 . 20 ) 


where {p{x ), a{x ), (5(a;)} are solutions of (2.5) with the identification 


cos X = y . 


( 2 . 21 ) 


3 Stability of level-0 boson stars 

We now proceed with the stability analysis of the boson star solutions (2.20) with 
respect to the linearized fluctuations partially breaking the R-symmetry to SO{5). To 



this end we assume, to linear order in A: 


0:3 = AA3(a;) dY^{6) cos(/ct), 0:4 = XA^i^x) dY^{6) sin(A;t), 


0:5 — 


sin^ xe 

cos® X 


1 + XA5{x) Yi{ 6 ) cos(/ct) 

g-5(x) 


9tx = Xf{x) Ye{6) sm{kt ), ci = 


a[x) 


cosx 


{l + \fi{x) Yi{ 9 ) cos{kt)) , 


C2 = —F==-(1 + A/2(a;) Yi{e) cos(/ct)) , C3 = tan a; (1 + Xfsix) Y^{e) cos{kt)) , 

\/a[x) cosx 

C4 = 1 + A/4(x) Yi{ 6 ) cos{kt) , C5 = 1 + A/5(x) Yi( 9 ) cos{kt ), 

0 = 01 + i02 = p{x) + \hi{x) Ygi^O) cos(/ct)^ exp [i (cat + A/i2(x) 50 ( 6 *) sin(A;t))] , 

(3.1) 

where 50 are the S'®-spherical harmonics, 


^5^50 = —s 50 — 


' + A)Ye, dY, = deY,. 


To order 0 {X), the 5 -form equations ( 2 . 19 ) are solved with 

^3 = ^4 = ^5 = 0, 04 = 05 = 0, 03 = “^01 “ ^02 • 


( 3 . 2 ) 


( 3 . 3 ) 
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Next, substituting (3.1) (with (3.3) and using the radial coordinate y, see (2.21) ) into 
(2.18) we hnd at 0{\): 


1 


4:ya?{y — l)p'k 
IQy^a^il - yY{p’Y + Sa^y‘^{l - yf{p^e^^uj‘^y{y - 1 ) 


A^p'a^py‘^{l — yY 

+ 2a(5a + y — 2)){p'Y + 4a|/(—1 + y){—2auYp^ — u‘^p‘^y + 2u‘^p‘^ + 

Yp‘^y‘^(l — yYe^^ — Aa^{aY{y — 1) + 4at'(|/ — 1) + + Gay + y"^ — 8a — Ay + A)') f 2 

1 


u 


/ 


48(1 - y)3l2y3i2a3pp>k I 4a y (1 - I/) (6p e u y - ai - Aai){p') + ay{i upy 

— Yoj'^p^ + Aioj'^p'^y — 2Aauj‘^p^ — Aiu'^p'^ — 12oj‘^p^y + 24a;^p^ + Gak^e^^ 

+ GpVu\l - y)e^^ + 2aHY + A){2a + y- 2)^ / + 

- 4|/V^(1 - yf{p'f + Aa^y{l - y){p'f + (pV|/(|/ - 1)6^'^ - 2pa{a -y + 2))p' ) hi 

' -lGy^a^{l - yfip'f + {8a^y^{y - 1) 

ph2e^^kuj 


A8p'a‘^py'^{l — yY 

+ 8aV(l - yf{a + 2y- A)){p'Y ) fi - 


Aya‘^{y — l)p' ’ 


(3.4) 


0 = h" + 


2ap'y{y — 1) + ap{2y — 1) — py + 2p , e^^y{Gp‘^u‘^ + k^ — «£(£ + 4) 

a{y - l)yp ^ A{y - l)a^y^ ^ 

-4|/V(i - yf{p'f + ^a^y{y - Yp' + p^^‘^y{y - i)e^^ 

hie^^ku 


ue 


25 


Ay^/‘^pa^{l — yY^‘^ 

+ 2pa{a + y-2)]f + 


2aYi-y)y' 


P r, (2a|/ -a-2y + A) 3h2Ujp^ + fik 
^ J ' -i\ J ' 


‘^ayiy - 1) 


(1 - yy/^yy^a ’ 


3p2a; , e‘^^y{Gp‘^uj‘^ + ky - ai{i + A) 1 

^ k ^ A{1 — yy/'^y^/'^ak 2a?{y — l)y 


+ <^Vy(y - l)e“ + 2a(o + 9 - 2) 1/2, 


(3.5) 

(3.6) 

-4(p')V(l - V)V 


(3.7) 
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n_ // , ^Vl/(l/-l)e^'^ + 2a(2a + |/-2)-4(p')^aV(l-|/)% , 

U — /2 + ^ _9/,. IN,. + 


Aa?{y - l)y 

uj‘^p^y{y — l)e^^ — 2a(3a — y + 2) — 4(p')‘^a^y‘^(l — y"^) 
Aa?{y - l)y 


4(1 — yy/'^c 


/a - 3pp'hi. 


(3.8) 


Additionally, there are two second order equations, which are consistent with (3.4)- 

(3.8), 

n X// , 4(p')V|/2(l-|/)2_a;yy(^_l)e2^-2a(a + |/-2)^, ^ 2ay-a-y + 2^., 

2 am-y)v <av(.V-l) 

^2Skyi/2f> ^25y^Qp2^2 _ p,2^ _ + 4) „ 3p^e^^u‘^y + 4a 3e^^p‘^kuh2 

~ 2{l-yy/^a 4{y - l)a^y^ ^ 2 (|/- l)aV + 2a?{y-l)y 

e^^k{4{p'ya^y‘^{l — yY — u‘^p^y{y — l)e^'^ — 2a{2ay — y + 2))f 3e^^hiu‘^p^ 

8yi/2(i _ |/)3/2 q3 ^2a^{y — l)y'’ 

(3.9) 


„ e^^ky^/^f 4{p'ya^y^{l - yf - u^p^y{y - 1)6^^ - 2a(2a + y - 2) , 

U 72+ , „.M/2„ r ,„ 0 /. 


2{l-yy/h 


2ayy-l)y 


, , ke^y4{p'ya^y^{l-yy-uj^p'^y{y-l)e^^-2a{2ay-3a-y + 2))f 

+- Symi - VY''V - 

2ay — y + 2j,, ^ e^^k'^y — a(? — 4ai — Sa ^ ,, 2 . 


ay{y - 1) 


-f2 + 


4{y - l)o?y‘^ 


72 - 6(7) hi. 


Further introducing 

hi = ky‘/^Hu h^ = u:y"^m, / = yry;j,<'+w2f., 

h = \y^‘*'"^ F,, /, = 1 


(3.10) 


(3.11) 


the spectrum of S'0(5)-invariant fluctuations about AO(6)-symmetric boson stars is 
determined solving (3.4)-(3.8) subject to the asymptotic expansions: 
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asymptotically at the AdS boundary, ie., ?/ —)■ 0+: 
- 3f - eu^ -8i + 3A;2 


H,=l 

H2 =Hl, + 
1 


12(« + 3) V + 0{V), 

3H^ ^e + - QP 

12(£ + 3) 


y + 0 {y ^), 


. _ _ _ Fl{e + U^-^k^ + U-2k^) 

^1- 2^0^ y + U{y), 


F, = - 


F^{i + e) 


8a2F^f + Fo¥ + F^ee - mp^k^ + 8a2F^l 


2 7^2 


2(£ + l) 8(£ + 3)(£ + l) 

+ IQF^e + 7Fo^£A;2 - F^k^ - 96plk^ + 24Fo^£ + 22Fl;e ] y + 0{y^), 


F=F^ + 


{e + 8£ - A;2 + 8)Fl 
4(£+l) 


y + Oiy"^) 


(3.12) 


at the origin of AdS, i.e., z = 1 — y ^ 0+; 

=F.« + + Kj) * 

+ 0(z^), 


(3.13) 


7/2 =Hi,„ + + Hi„i 

n =do + + ^f;:/ 

+ f;;„«+ jdo)‘ + o(^"). 


(3.14) 


(3.15) 


+ 4do(<iS)"(p5)"‘^" + - ^F,%f - ^F;‘P - iffoj ^ + O(z^) , 


(3.16) 
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5 


F =^Koip,y^ + 2^10 + i^^cj^(psrff2y +w - gff2ypyyd^r 
- §fU4)W^^ - :^(4 )Vf!‘_„ + + 3hIM)V 


24 


24 


12 


- - ^(4)"H2Vp5y^ + o(^') ■ 

(3.17) 

Note that without loss of generality we normalized linearized fluctuations so that 
Hi =1. Furthermore, the total number of parameters characterizing the solu- 

y=0 

tion, 

{P, , F >, , F,\} , (3,18) 


is precisely the order of the ODE system (3.4)-(3.8). 

Once again, we use the shooting method developed in [14] to determine the spectrum 
of the linearized fluctuations. First we outline the computations for the light boson 
stars, i.e., for po 1, and then present the results for the general boson stars. We 
restrict out attention to j = 0 (ground state) boson stars. For j > 0 the radial prohle 
of a boson star p(y) has j nodes inside the interval y G (0,1). This results in additional 
poles in the fluctuation equations (3.4) and (3.5) which renders our shooting method 
inapplicable^. 


3.1 Spectrum of linearized fluctuations of j = 0 light boson stars 

Using (2.11), to leading order in pq, 

piy) = Pj=o{,y) = Po + 0 {pI) , u = ujj=o = A +0{pl), (3.19) 

and 

H, = + 0 {pI) , + 0 {pI) , 

F = py^iy) + 0{pt ), Fi = plF,^ 2 {y) + 0{pI) , F^ = pIF^M + 0{pl ), 

(3.20) 

we hnd from (3.4)-(3.8): 

^ _ tt! ^yFl^ 2 ,o . + F) — k'^y + y + 3)F2^2 (^^ — 1)-^ i ,2 

'’°+F(|/-l)+ 2 Ay^{y-l)k^ ^ 2 A {1 - y)y^k^ 

{ 2 i‘^y + 8 iy + 6 k‘^y)F ^2 iijy - 1 ) - y - l)hfi,o i 2 iy - 2 k‘^y)H 2 fl 
96?/3/(;2(l - y) 2 y{y - 1) yk‘^{y - 1) 

^Identical technical difficulties were also observed in [1]. These difficulties can be resolved assuming 
a more general fluctuation ansatz for j > 0 as in [15]. 
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^ , (£(|/- 1 ) + 5 |/- , {e + U-k^)H2,, 

U =^2,0 -^-^ 


y{y-^) 


My-^) 


My - 1) ’ 


0 =f;, + 


{iy- i + MM M 2 4:8yH2,o 


My-'^) 


y{y-^) y-i 


(3.22) 

(3.23) 


0 —2 + 48 ?/ H2 Q + 


2„, , {e-ey + Ai)F2 , (^ + 4)Fi,2 


Ay 


+ 


2// 


+ F2,2 + 24?/i72,o^ (3.24) 


n T?> , M2 , ( 2 £(?/- 1 ) + 10?/-2)F2,2 , k‘^F^2 ,2 

0 =^2 2 + 7T H--7^-h — ^- 6iii nfc 


(3.25) 


2y{y - 1) ' Ay{y - 1) 'Ay 

The system of ODEs (3.21)-(3.25) can further be reduced to a single 4th-order ODE 
for 

, , 2(£?/-£ + 7?/-4)^,„ , 1 

U = 4/1,0 --D- 4 (i,o + 


yiy - 1) 


2 //^(i - yf 


3iM - 5^"?/ + 3Qiy^ - k'^y^ + 2i'- 


2„.2 


- 50iy + k^y + 9Ay^ + 14£ - 106?/ + 24j H'l, + 

- 13 f + M + 60iy - bk^y - AOi + 3^ + 50y - 30^ + ^Qy 2 (l _ yy + 

- 2fk^ + 6Af - me + e- QAkA . 
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Solving (3.26) with the boundary conditions 


(3.26) 


(3.27) 


(3.28) 


lim Hifi{y) = 1 > Mfiiy) = fiiiite , 

y^0+ 

determines the perturbative spectrum e. 

Note that (3.26) with (3.27) allows for a polynomial-in-?/ solution: 

HiM) = 4^i,oy,n,±(//) = 24^1 {-n ,/ + r? + 4;/ + 3;?/), 
e = kl^ _ = (£ + 2nf , e = kl^^^ = (£ + 8 + 2??)^ , 

where n = 0,1, • • ■ indexes the excitation level of a (light) ground state boson star; 
additionally ± denotes discrete branches for a hxed n. Although for different branches 
44i,o,£,n,- = 44i,oy,n,+, the radial prohles for the remaining fluctuations do differ; e.(/., 
for n = 1, solving (3.22) with the appropriate boundary conditions, (3.12) and (3.14), 
we End for 772 ,0 = 472,0,?:,n,±, 
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(3.29) 
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Figure 2: Low-lying states in the fluctuation spectrum of j = 0 boson stars as 
a function of charge Q: £ = {0,1} (left panel, right panel), n = {0,1,2,3} 
(blue,red,green,magenta). Leading order in Q —)■ 0 frequency eigenvalues kfn- are 
given by (3.28). 

3.2 Fluctuation spectrum of j = 0 boson stars 

Perturbative in pq spectrum of the j = 0 boson stars (3.28) provides a starting point for 
the construction of the generic states labeled by {£, n, ±}. These states are obtained 
numerically solving (3.4)-(3.8), subject to the asymptotic expansions (3.12)-(3.17). 
The results of the analysis are presented in figure 2. States with £ = 0 represent 
S'0(6)-invariant fluctuations. Furthermore, the state {£ = 0,n = 0, —} is a zero mode 
corresponding to rescaling of A in (3.1). We observe that the state {£ = l,n = 0, —} 
becomes tachyonic for large values of the boson star charge Q (blue line, right panel). 
For small values of the charge Q the frequency eigenvalues of the fluctuations are close 
to kjn^± (given by (3.28)) and thus are stable. 

4 Conclusion 

The low-energy dynamics of strongly coupled gauge theories in a finite volume is rather 
involved. Motivated by the gauge theory/string theory correspondence we discussed 
stability of stationary solutions, boson stars, supported by a complex scalar held in 
AdS^ X S^. These solutions are typical representatives of S'0(4)-invariant states in 
global AdS^ that fail to gravitationally collapse in the limit of vanishing mass [1,2, 
5]. Unlike small smeared black holes in AdS^, which are unstable with respect to 
localization on [8,9,11], we explicitly demonstrated that boson stars are stable with 
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respect to linearized perturbations leading to spontaneous breaking of the global SO{6) 
symmetry below some critical mass (an in particular in the limit of vanishing mass). 
The result is far from being unexpected: the Gregory-Laflamme instability is triggered 
by the localization of energy in a small region of the space-time, smeared over a large 
compact transverse space. A characteristic feature of a boson star (and also generic 
solutions that fail to gravitationally collapse in global AdSd+i) is a broad distribution 
of the matter profile — quite an opposite regime for the GL instability. While the GL 
instability excludes low-energy R-symmetry singlet equilibrium states in A/” = 4 SYM, 
our analysis suggests that there is a large class of low-energy intrinsically nonthermal 
states in the theory, invariant under the R-symmetry. 

Our discussion was done in a toy model of the gauge theory/gravity correspondence. 
It would be interesting to study boson stars in type IIB supergravity proper. A starting 
point in this direction could be the supergravity solutions used in models of holographic 
superconductors, as in [16]. 
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